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Abstract. We present sets of equation of state (EOS) of nuclear matter including 
hypcrons using an SUy(3) extended relativistic mean field (RMF) model with a wide 
coverage of density, temperature, and charge fraction for numerical simulations of core 
collapse supernovae. Coupling constants of E and S hyperons with the a meson are 
determined to fit the hyperon potential depths in nuclear matter, Us(po) ~ +30MeV 
and Us(po) — — 15MeV, which are suggested from recent analyses of hyperon 
production reactions. At low densities, the EOS of uniform matter is connected 
with the EOS by Shen et al., in which formation of finite nuclei is included in the 
Thomas-Fermi approximation. In the present EOS, the maximum mass of neutron 
stars decreases from 2.17M Q (Ae/z) to 1.63M Q (NYefi) when hyperons are included. In 
a spherical, adiabatic collapse of a 15Af© star by the hydrodynamics without neutrino 
transfer, hyperon effects are found to be small, since the temperature and density do 
not reach the region of hyperon mixture, where the hyperon fraction is above 1 % 
(T > 40MeV or p B > 0.4 fm~ 3 ). 



1. Introduction 

The equation of state (EOS) plays an important role in high density phenomena such 
as high energy heavy-ion collisions, neutron stars, supernova explosions, and black 
hole formations [DI2ll3lllll5lEllZllSllSlinainilI2llI3llI3. The recent discovery 
of the strongly interacting quark gluon plasma (sQGP) [14] attracts attentions to the 
EOS and transport coefficients in the quark gluon plasma (QGP). The core region 
of neutron stars, where matter becomes very dense (~ 10 15 g/cm 3 ), is an interesting 
play ground of quark and hadronic matter models. Various ideas for the new form 
inside neutron stars have been proposed including strangeness and quark degrees of 
freedom [TJ, [2J, [3J, SI El El El El E] - Core-collapse supernovae also involve high density 
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and temperature. The nuclear repulsion at high densities drives the shock wave at core 
bounce and the passage of shock wave heats up the matter inside the supernova core. 
A hot, lepton-rich neutron star (proto-neutron star) is born after the explosion and 
cools down by emitting supernova neutrinos. When black holes are formed for more 
massive cores, extremely high density and temperature are involved, where hyperons 
should appear and quarks would be deconfmed. In order to describe the whole evolution 
of core-collapse supernovae by numerical simulations, one needs to prepare the set of 
microphysics under such extreme conditions. One of the most important ingredients is 
the set of equation of state (EOS) that contains necessary physical quantities. It is to 
be noted that one must cover a wide range of temperature, density and composition in 
a consistent manner and theoretical framework. 

Until now, the two sets of EOS (Lattimer-Swesty EOS [15] and Shen EOS [16]) have 
been widely used and applied to numerical simulations of core-collapse supernovae pUl 
ITT] and black hole formations [12], [13]. The Lattimer-Swesty EOS is based on a 
compressible liquid-drop model, whose mass and mean field potential are motivated 
by non-relativistic zero-range Skyrme type interactions. The Shen EOS is based on a 
relativistic mean field (RMF) model, whose interactions are determined by fitting the 
binding energies and nuclear radii of stable as well as unstable nuclei [17]. Coexistence 
of nuclei and uniform matter is included in the Thomas-Fermi approximation in 
the Wigner-Seitz cell, and the alpha particles are assumed to follow the statistical 
distribution with excluded volume effects. 

The constituents in these EOSs are neutrons, protons, alpha-particles and nuclei, 
restricting the framework within the non-strange baryons. These degrees of freedom 
may be enough to simulate the early stage of hydrodynamical evolution of supernova 
explosions. However, in order to clarify the long-time evolution from core collapse [11] 
to proto-neutron star cooling [18], [19], black hole formation [121 US] , neutron star 
mergers and gamma ray bursts that may involve higher density/temperature, it would 
be necessary to include other particle degrees of freedom. Especially, hyperons (baryons 
containing strange quarks) are commonly believed to appear in neutron star core and 
to modify the neutron star profile [H [21 [31 HI [5], El [7J. While there are several works 
which include the hyperons in the proto-neutron star cooling [19], there has been no 
study on the dynamics of core-collapse supernovae adopting the EOS with hyperons. 
This is partially because EOS table of supernova matter including hyperons has not 
been available in public. In addition, the determination of the interaction for hyperons 
has been difficult having large uncertainties so far. 

Recently, developments in hypernuclear physics have narrowed down the allowed 
range of hyperon potential depth in nuclear matter. The potential depth of A has 
been well known to be around (po) — — 30MeV from bound state spectroscopy. 
For E baryons, it was considered to feel similar potential to A, because it contains 
the same number of light (u, d) quarks. From the recently observed quasi-free £ 
production spectra [20] , it is now believed that E baryons would feel repulsive potential 
in nuclear matter; £7x(p ) ~ +30MeV J2TJ [221 122], Also for E baryons, the analyses of 
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the twin hypernuclear formation [21] and the S production spectra [221 [2SJ [21] , suggest 
the potential depth of around U~ N \po) ~ — 15MeV. These £ and 5 hyperons are 
particularly important in neutron stars, since nuclear matter can take a large energy 
gain from neutron Fermi energy and symmetry energy by replacing, for example, two 
neutrons with a proton and a negatively charged hyperon (E~ or The updates 

on the interactions of hyperons may have impact on supernova dynamics and thermal 
evolution of proto-neutron stars. 

In this paper, we present new sets of EOS of dense matter with hyperons, 
abbreviated as EOSF, under the current understanding of interaction. We provide the 
data table covering a wide range of temperature (T), density (p B ), and charge-to-baryon 
number ratio in hadronic part ( Yq) , which enables one to apply to supernova simulations. 
Our framework is based on the RMF theory with the parameter set TM1 [17], which was 
used to derive the EOS table by Shen et al. [16], and is extended to include hyperons 
by considering the flavor SU(3) Lagrangian [2j. Therefore, our EOS table is smoothly 
connected with Shen EOS and can be used easily as an extension of Shen EOS table in 
numerical simulations. 

It is well known that the RMF predicts large values of incompressibility (K ~ 
300 MeV) and symmetry energy, and these are sometimes considered to cause problems 
in applying to dense matter EOS, since they lead to too high maximum mass of neutron 
stars without hyperons and may be unfavorable to core-collapse explosions. It should 
be noted that K and symmetry energy values are not yet well determined separately in 
a model independent manner. Analyses of collective flow data at AGS energies suggest 
K = 210 — 300 MeV [281 129], and collective flows at SPS energies are shown to be 
more sensitive to the mean field of resonance hadrons rather than to the cold matter 
EOS [30]. For symmetry energy, it is possible to describe binding energies, proton- 
neutron radius differences and isovector giant monopole resonances simultaneously by 
incorporating density dependent coupling (DD-ME1) in RMF and relativistic RPA [31] 
with a larger value of symmetry energy than that in non-relativistic models. Note 
that the interaction of RMF-TM1 is constrained by the nuclear masses, radii, neutron 
skins and excitations [XT] . The large value of K leads to a large neutron star mass 
(stiff EOS), which seems not preferable for explosion. On the other hand, the large 
symmetry energy is known to be preferable for explosion, having less free proton 
fraction and smaller electron captures. These two effects are competing each other 
in the sophisticated numerical simulations pj]. We note also that the RMF fulfills 
automatically the causality (the sound velocity should not exceed the light velocity) 
whereas the non-relativistic frameworks breaks down at high densities appearing in the 
simulations. Thus at present we do not find problems in applying RMF EOS to dense 
matter compared to non-relativistic models. 

This paper is arranged as follows. In section [2} we describe the framework to 
calculate the dense matter at finite temperature including hyperons. We explain the 
updated information on hyperon potentials in nuclear matter and adopted potential 
values in EOSF. We describe also the prescriptions to provide the data table for 
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the wide range of density including sub-saturation densities where finite nuclei appear. 
In section [3j we report the properties of EOSF in comparisons with nucleonic EOS 
(TMl/Shen EOS). We apply EOSF to cold neutron star matter and supernova matter. 
We show several properties of EOSF at finite temperatures by examining energies, 
chemical potentials and compositions. The data tables are successfully applied to 
hydrodynamical calculations of adiabatic collapse of iron core of massive stars. We 
examine the possibilities of hyperon appearance in supernova cores. Summary and 
discussions are given in section HI In the appendix, we provide the the definitions of 
quantities in EOSF. 

2. Model and Method 

In this work, we construct the EOS table of supernova matter based on a relativistic 
mean field model. We adopt the parameter set TM1 [17] for non-strange sector. For 
its flavor SU(3) extension, we start from the work by Schaffner and Mishustin [2], and 
we include the updated information on hyperon potentials from recent experimental 
and theoretical hypernuclear physics developments. Low density part of the EOS is 
connected with the Shen EOS. 

2.1. Relativistic mean field model with hyperons 

The relativistic mean field (RMF) theory is constructed to describe nuclear matter and 
nuclei based on the relativistic Briickner-Hartree-Fock theory [32], which successfully 
describes the nuclear matter saturation. It is preferable to adopt the relativistic 
frameworks for astrophysical applications, since they automatically satisfy the causality, 
i.e. the sound velocity is always less than the speed of light. 

The RMF parameter set TM1 is determined to describe binding energies and 
nuclear radii of finite nuclei from Ca to Pb isotopes and fulfills the nuclear matter 
saturation. The incompressibility of symmetric uniform matter and the symmetry 
energy parameters are found to be K = 281 MeV and a sym = 36.9 MeV. When it 
is applied to neutron stars, the maximum mass of cold neutron stars with TM1 is 
2.17M . 

The extension of the RMF to flavor SU(3) has been investigated by many authors. 
A typical form of the Lagrangian density including hyperons is given as [2], 

£ = (ip - M B ) m B + \d^d^ - U a (a) 

B 

- -o^uv + ^ml^u, - ■ R» v + l -m 2 p & ■ R, 

B 

U a (a) = -m a a + —a 6 + —a , 
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c yy = \d v Qd v c - l -mie - \^r u + 

- ^ (g (B C + m^v) , (l) 

B 

where the sum runs over all the octet baryons. In this Lagrangian, hidden strangeness 
(ss) scalar and vector mesons, ( and <fr, are included in addition to a, uj and p 
(represented by mesons. Strength tensors of u, p and (ft mesons are shown in u^, 
R^ u and respectively. The Lagrangian contains meson masses, coupling constants, 
and self-coupling constants as parameters. 

In introducing hyperons in RMF, we have large ambiguities in hyperon-meson 
coupling constants. One of the ways to determine the parameters is to rely on 
symmetries. Schaffner and Mishustin [2] have determined hyperon-vector meson 
coupling constants based on the SU(6) (flavor-spin) symmetry, 
111 1 

gfiW — — 2 guS ~ ^ uB ' 9pN ~ 2^ ps ~~ 9pE ' 9 p a — , (2) 

2g$A = 2^ = g</>E = — — 9ojn , g<t>N = . (3) 

Scalar mesons in RMF may partially represent contributions from some other 
components than qq, such as irn in a. In Ref . [2] , the scalar meson couplings to hyperons 
have been given based on the assumption that hyperons feel potentials in nuclear and 
hyperon matter as, 

U[ N) = U^ N) = -30 MeV , ui N) = - 28MeV , 

U P ~ u a ] ~ U P ~ 2U a ] ~ -4° MeV , 
where U B denotes the potential of B in baryonic matter at around po composed of B ' . 
Recent developments in hypernuclear physics suggest that hyperon potentials in nuclear 
matter are repulsive for S [201 1211 [221 123] , and weakly attractive for 3 1251 126| 127] . 
respectively. 

3 hyperons are expected to have nuclear bound states, and the bound state 
spectroscopy at forthcoming facilities such as J-PARC and FAIR will give a strong 
constraints on the 3 potential in nuclear matter. At present, the depth of the 3~-nucleus 
potential has been suggested to be around 15 MeV from the analysis of twin hypernuclear 
formation [24] and the (K~ , K + ) spectrum in the bound state region [25j. In the former, 
the binding energy of the 3~-nuclear system is found to be consistent with a shallow 
3~-nuclear potential in an event accompanied by two single hyperfragments emitted 
from a 3~ nuclear capture at rest (a twin hypernuclei) found in a nuclear emulsion |24j . 
In the latter, while the resolution of experimental data is not enough to distinguish the 
bound state peaks, the observed yield or the spectrum shape in the bound state region is 
found to be in agreement with the calculated results with ~ —15 MeV [25l |26| [27] . 

For E hyperons, it is necessary to analyze continuum spectra. In the observed 
(quasi-)bound £ nucleus |He [33], the coupling effect is strong and the repulsive 
contribution in the T = 3/2, 3 Si channel is suppressed, then it does not strongly 
constrain the S potential in nuclear matter. The analysis of S~ atomic data suggested a 
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E~-nucleus potential having a shallow attractive pocket around the nuclear surface and 
repulsion inside the nucleus [341 [35], [36] , but the atomic energy shift is not sensitive to 
the potential inside the nucleus. In the distorted wave impulse approximation (DWIA) 
analyses of the quasi free (QF) spectrum in the continuum region [20j (2TJ [221 123], it is 
suggested that the E hyperon would feel repulsive real potential of 10 ~ 90MeV. Recent 
theoretical analyses favor the strength of repulsion of around +30MeV [211 [22], [23] . This 
repulsion may come from the Pauli blocking effects between quarks due to the isovector 
nature of the diquark pair in E [37]. In a Quark-Meson Coupling (QMC) model, medium 
modification of the color hyperfine interaction in the quark bag is found to be the origin 
of repulsive E potential [38]. The E potential in nuclear matter at saturation density 
is predicted to be around +30 MeV (repulsion) in a quark cluster model YN potential 
[37], and a chiral model also predicts a similar repulsion [39] . 

From these discussions, we adopt the following potential strength as recommended 
values, 

U^ N) (p ) ~ +30 MeV , ui N) (p ) ~ -15 MeV . (4) 

The above spectroscopic studies have been done mainly with non-relativistic 
frameworks for hyperons, then the potential should be regarded as the Schrodinger 
equivalent potential in RMF. The Schrodinger equivalent potential is related to the 
scalar (U s ) and vector (U v ) potentials as, 

U B ( P ,E(p)) = U s (p) + ^U v (p) 

= 9*bo- + g (B ( + -tj (uw + 9 p bR + g<f,B<t>) , (5) 

where R represents the expectation value of the p meson. We have fixed g a B value by 
fitting the hyperon potential depth in normal symmetric nuclear matter, 

ui N \po) = 9«BO- {N \p Q ) + g. B ^ N \po) , (6) 

where a^ N \po) and u^ N \po) represent the expectation values of a and to mesons in 
symmetric nuclear matter at po- We adopt the parameter set TM1 for nucleon sector, 
and we determine g a ^ and g a z to reproduce the potential depths of E and H hyperons in 
Eq. (T4]) as listed in Table [U We choose other hyperon-meson coupling constants referring 
to the values in Ref. [2]. We show the values of g a ^, and g a s for different potentials in 
Table [SJ 

In the literatiure, the instatibility due to the negative effective mass of nucleon 
has been reported [4*01 [41] . As o increases, the nucleon effective mass reaches zero at 
a = Mjv/^o-jv where hyperon effective masses are still positive and act to further increase 
a, leading to the nucleon negative effective mass. This effect depends very much on aY 
couplings, which are small within the current sets of parameters, and we did not find 
this instability in the range of data tables. However, we found this instability occurs at 
very high densities/temperatures which are relevant in the black hole formation [42]. 
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m a (MeV) g 3 (MeV) g 4 CuJ 
511.198 1426.466 0.6183 71.3075 



9mb 


a 


c 




P 





N 


10.0289 





12.6139 


4.6322 





A 


6.21 


6.67 


8.41 





-5.95 


£ 


4.36 


6.67 


8.41 


9.26 


-5.95 




3.11 


12.35 


4.20 


4.63 


-11.89 



Table 2. The coupling constants of SiV and ^N. 



U^\p ) (MeV) 


9<rE 




+90 


2.58 




+30 


4.36 


present 





5.35 




-10 


5.63 




-30 


6.21 


Ref. [2] 


t/f j (p ) (MeV) 


g*~ 




-15 


3.11 


present 


-28 


3.49 


Ref. [2] 



2.2. Free thermal pions 

In black hole formation processes as found in Ref. [12] , the temperature goes up to 
around 100 MeV. At these high temperatures, pion contributions become dominant. 
Charged pions may condensate at high densities in neutron star matter [HE]- ^° 
estimate the effect of pion mixtures, we also prepare the EOS table including free 
thermal pions assuming the pion mass is not affected by the interaction. This is of 
course oversimplification, however, the first trial to include pions. Further sophisticated 
studies are necessary. 

The density of free thermal pions is calculated to be 



(2tt)3 exp((^(p) - fi^/T) - 1 ' w 

where p n = pc, 0, — pc for 7r + ,7r°,7r~, respectively. When the absolute value of the 
chemical potential reaches the pion mass, pion condensation occurs; i. e. the amount of 
condensed n at zero momentum can take any value at pc — We have determined 

the amount of condensed tt in the following way. First we solve the equilibrium condition 
and obtain p B and pc without condensed ir. When \p c \ > m^, we set \pc\ = m n 
and re-evaluate hadron densities except for the condensed n to satisfy the condition of 



= p°:° nd + 
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Pb = Pb (Given). Finally, the amount of condensed 7r is given so as to satisfy the charge 
density condition, p c = Y c ( Given) p B . 

The pion condensation in the current treatment is a simple s-wave Bose-Einstein 
condensation, which is different from the pion condensation derived from p-wave nN 
interaction [8] . We mention that pion condensation will be suppressed after considering 
s-wave nN repulsive interaction discussed in the energy shift of deeply bound pionic 
atoms [131 HH H5] and in pion-nucleus scattering [46J . 

2.3. Low density 

By using these potentials, we can immediately obtain a EOS of uniform dense matter 
with strangeness based on the RMF theory. We also need to cover the low-density region 
below po, where the inhomogeneous matter appears. Here we connect the uniform matter 
EOS with Shen EOS [16], which is based on the same RMF parameter set TM1 and 
treats the inhomogeneity with the Thomas-Fermi approach. 

We include the contribution from inhomogeneity by adding the free energy 
difference of Shen EOS values from those in uniform matter, 

F = F][ MF + AF Nud (8) 

at p B < po, where 

AF Nuc i = F S hen — Frmf ■ (9) 

Other variables are derived in the same way as the above equations. The deviation due 
to inhomogeneity AF vanishes at p B > po- These prescriptions produce the extended 
EOS tables for studies in astrophysics, containing inhomogeneity at low density and 
strangeness information. The compositions of n,p, a, A, Y are consistent with Shen 
EOS table and the sum of each component ratio becomes unity. 

2.4- Tabulation of thermodynamical quantities 

Thermodynamical quantities are provided in the data table as a function of baryon mass 
density p B , charge ratio Yq, and temperature T. Here Yq means the charge ratio defined 
as Yc = nc /tib and ric is a charge density. See Appendix for the list of quantities and 
their definitions, which are slightly revised from the original table of Shen EOS. 

For the purpose of numerical simulations, we prepare the EOS table containing the 
contributions of leptons and photons by adding the energy, pressure and entropy from 
electrons, positrons and photons to the hadronic EOS. We treat electrons and positrons 
as ideal Fermi gas with the finite rest mass and calculate photons according to the 
standard expressions for radiations. 

The baryon mass density, charge ratio and temperature cover the following range, 

• p B = 10 5 - 1 ~ 10 15 - 4 (g/cm 3 ) (104 points) 

• Y c = and 0.01 ~ 0.56 (72 points) 

• T = and 0.1 ~ 100 (MeV) (32 points) 
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Mesh points for p B , Y c (> 0) and T(> 0) are taken as approximate geometric sequences 
with Alog 10 p B = 0.1, Alog 10 Y c = 0.025 (-2.00 < log 10 V c < -0.25) and Alog 10 T ~ 
0.1, respectively. These ranges and mesh points are the same as those in Shen EOS. By 
connecting smoothly in the way described above, the EOS table is combined with Shen 
EOS at lower densities below p Q while it includes full baryon octet at high densities 
so that one can see the effects of hyperon mixture. Some tabulated quantities in the 
EOS table need attentions. The values of the mass A and charge number Z of heavy 
nucleus are taken from Shen EOS at densities below po and are set to be zero above po- 
Similarly, the fraction of a-particle and heavy nucleus are taken from Shen EOS at low 
densities and are set to be zero at high densities. 

3. Properties of EOS tables and astrophysical applications 

We report the properties of dense matter in the present EOS table with hyperons 
(EOSY) and their applications to neutron stars and supernovae. We adopt hereafter 
the case of (Uj- \(po),Ua ■ (po)) = (+30 MeV, — 15 MeV) as a standard case, which 
is currently the most recommended set of hyperon potentials. We also consider the 
case with pion contribution (EOSVtt) and the attractive hyperon potential case [2] 
(U^ f) (p Q ),U^ f) (p Q )) = (-30 MeV, -28 MeV), abbreviated as EOSV(SM). 

3.1. Neutron star matter 

We first study the EOS of neutron star matter, which is under the (3 equilibrium at zero 
temperature. We here add electron and muon contributions under the (3 equilibrium 
and charge neutrality conditions. We consider uniform matter ignoring finite nuclear 
effects. 

We show particle compositions in neutron star matter in Fig. [T]to see the appearance 
of new degrees of freedom. We display the cases of nucleonic (TM1, upper-left) and 
hyperonic (EOSV, upper-right) EOS. Results with hyperonic EOS with attractive E 
potential (EOSV(SM), lower-left) and hyperonic EOS with pions (EOSV7T, lower-right) 
are also shown for comparison. The particle composition of neutron star matter is very 
sensitive to the choice of hyperon potentials. With attractive E potential, E~ appears 
at lower densities than A. With repulsive E potential, A appears first followed by 
H~ and S°. This behavior is different from the previous works that adopt attractive 
potentials P EJ El 13 E], and pointed out in Refs. [21 EJ HJ EH]. When we allow 
the appearance of pions, condensed pions (ir~) appear prior to hyperons. With n c 
condensation, the charge chemical potential is restricted to be |p c | < m n and the proton 
fraction becomes larger, then the neutron chemical potential is reduced. As a result, 
the threshold density of hyperons are shifted up. The density of A appearance is about 
0.37 fm~ 3 and other hyperons such as S~ are also suppressed. 

In the left panel of Fig. [2l we show the energy per baryon (E/B) and chemical 
potentials (p n and p p ) in neutron star matter in nucleonic EOS (TM1), EOSV and 
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Table 3. Constituents, assumed hyperon potentials, threshold densities, maximum 
masses of neutron stars, central densities giving maximum masses of neutron stars and 
neutron star masses at the threshold central densities in nuclconic EOS (TMl/Shen 
EOS [HI US]), hyperonic EOS with attractive hyperon potentials (EOSF(SM) [5]), 
hyperonic EOS with repulsive hyperon potentials (EOSK, present work) and hyperonic 
EOS with repulsive hyperon potentials including pions (EOSy-7r, present work). For 
7r~, the maximum density of condensation is also shown in parentheses. Threshold 
densities of protons and muons are 1.1 x 1CP 4 fm~ 3 and 0.11 fm -3 , respectively. 



EOS 


TMl/Shen EOS EOSF(SM) 


EOSF 


EOSF7T 


Constituents 


Ne(p) 


NYe(p) 


NYe(p) 


A^y7re(;u) 


J 7 (AM /Tl IT T T\ 

Uj, ' (MeV) 




-30 


+30 


+30 


Uk ' (MeV) 




-28 


-15 


-15 


p {thr) (fm -3 ) 










A 




0.32 


0.32 


0.37 


E" 




0.29 


1.14 


1.1 


S° 




0.57 


1.34 


1.3 


S+ 




0.69 


1.47 


1.5 






0.43 


0.40 


0.56 


wo 




0.62 


0.71 


0.74 










0.16(0.88) 




2.17 


1.55 


1.63 


1.65 


pL max) (fm" 3 ) 


1.12 


0.79 


0.79 


0.97 


M^\M & ) 




1.17(E-) 


1.28(A) 


1.22(A) 










0.51(71--) 



EOSYV. Compared with nucleonic EOS, E/B is much lower in EOSF at high densities. 
Chemical potentials are also suppressed with hyperons correspondingly. There are 
several origins for this energy gain. First, nucleon Fermi energy decreases due to 
the hyperon mixture. Secondly, the lepton contribution is suppressed when negatively 
charged hyperons emerge. In addition, the repulsive vector potential becomes small, 
because the ujY couplings are smaller than uN and the isospin asymmetry becomes 
smaller when negatively charged hyperons appear, as shown in the long-dashed lines in 
Fig. [3 

For E/B and /i n , pionic effects are small and only visible around p B ~ 0.4 fm -3 , 
while we find large differences in p p . The equality pc = Pp ~ Pn = ~ m n under 7r _ 
condensation reads the Fermi energy relation, Ep{n) = Ep(p) + m n . Since a neutron 
on the Fermi surface is replaced with a proton and a tt~ having the same total energy, 
we have to pay the cost of the pion rest mass energy in exchange for the Fermi energy 
reduction and symmetry energy gain. At higher densities where hyperons appear, pionic 
effects becomes smaller, and disappear at p B = 0.88 fm~ 3 . 

The s-wave pion condensation would be suppressed when we include the 7rA^ 
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Neutron Star Matter 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

Pb ( fm " 3 ) Pb ( fm " 3 ) 



Figure 1. Composition of neutron star matter in nucleonic EOS (TM1, upper-left), 
hyperonic EOS with attractive potential (EOSF(SM), lower-left), hyperonic EOS with 
repulsive potential without (EOSK, upper-right) and with pions (EOSy7r, lower-right). 
The number fraction of particles are plotted as functions of baryon density. The species 
of particles are denoted as in the legend. 



interaction. We evaluate the pion energy by using the potential of the form 

SHI HE], 



(b p B + hSp) + 1 + — =- Re£? p 2 B 



M N J v u ' ' ' V 2Mtv 

'free 



,(io) 



where b± = b^ f (l — ap B /po), 5p = p n —p P - As typical examples, we adopt the parameter 
sets from the analyses of pionic atom data; b = — 0.023/m^, b\ = — 0.085/m,,- (a = 0), 
ReB = -0.021/ (less repulsive) gl], and b = -0.0233/ttv, ^ rcc = -O.U73/m n , 
a = 0.367 {pi = — 0.1149 /m n at p B = 0.6po), HeB = — 0.019/m^ (more repulsive) 



In Fig. [31 we show the pion energy, E n = \/m\ + 2m 7T U s , calculated with these potentials 
and proton fraction in TM1 EOS. With these potentials, we find that the existence 
of the s-wave pion condensed region, E n < p e , depends on the pion optical potential 
parameters. Since the pion potential above the normal nuclear density is not yet known, 
the realization of pion condensation may be marginal and mo del- dependent. 

We apply the above four EOSs of neutron star matter discussed above (TM1 EOS, 
EOSy(SM), EOSy, EOSF7t) to the hydrostatic structure of neutron stars by solving 
the Tolman-Oppenheimer-Volkoff equation. We plot the gravitational mass of neutron 
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Figure 2. EOS of neutron star matter and supernova matter for nucleonic EOS 
(TMl/Shcn EOS), hyperon EOS(EOSY), hyperon EOS with free pions (EOSYtt). 
The upper panels display the energy per baryon, and the lower panels the chemical 
potentials of proton and neutron. The dashed line is for the nucleonic EOS, the solid 
line is for hyperon EOS (EOSY), and dotted line is for hyperon EOS with pions 
(EOSYtt). 
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Figure 3. Meson expectation values and electron chemical potential in neutron star 
matter for TM1 EOS (thin lines) and hyperon EOS(EOSY) (thick lines). Thin solid 
lines show the pion energy, E„ = y/m^ + 2m ir E/ s (7r _ ), where lower less repulsive (upper 
more repulsive) line shows the results with the pion optical potential in Ref. [H] ( |45j ) 
by using the proton fraction in TM1 EOS. 
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Figure 4. Neutron star masses are shown as functions of central density. The dashed, 
thin solid, thick solid, dotted lines show the results with nucleonic EOS (TM1 [17]). 
hyperonic EOS with attractive hyperon potentials (EOSF(SM) [2]), hyperonic EOS 
with repulsive hyperon potentials (EOSF) and hyperonic EOS with pions (EOSF7r), 
respectively. Filled points show the maximum neutron star masses, and open circles 
show the threshold densities of hyperons and pions. 

stars as a function of central baryon mass density in FigJH The maximum mass of 
neutron stars in EOSF is smaller than the case of nucleonic EOS because of the softness 
from hyperons. The maximum mass is 1.63 M & for EOSF in contrast to 2.17 M Q 
for nucleonic EOS when we adopt repulsive potentials for hyperons. The maximum 
mass is further reduced to be 1.55 M© in EOSF(SM) with attractive potentials. The 
neutron star masses with EOSFvr are reduced in the mid range of central density 
0.16 fm -3 < p B < 0.8 fm -3 , but the maximum mass (1.65M ) is almost the same. This 
is because the maximum mass is mainly determined by the EOS at densities around 
0.8 fm~ 3 or more, where the condensed pion density is small. The central density of 
a typical neutron star having 1.4M is 0.35 fm -3 in nucleonic EOS (TM1), which is a 
little above the threshold density of A in EOSF. In this case, hyperons are limited only 
in the core region, and the neutron star mass does not get a large reduction as seen 
in Fig. HI We summarize the neutron star masses and the threshold densities to have 
hyperons in neutron star matter in Table [3l 
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Figure 5. Composition of supernova matter at (T,Yc) = (10 MeV, 0.4) (left) and 
(T,Y C ) = (10 MeV, 0.2) (right) in the hyperonic EOS table without (EOSV, upper) 
and with (EOSy7r, lower) pions. The number fraction of particles are plotted as 
functions of baryon density. The species of particles are denoted as in the legend. 



3.2. Hyperonic matter at finite temperatures 

We next study the EOS of supernova matter, where the hadronic charge fraction (Yc) 
is fixed at finite temperature. We here show the results including electron and photon 
contributions. Finite nuclear effects in Shen EOS are included. The treatment of leptons 
in supernova matter are explained in Appendix A. 4. 

In order to demonstrate the contents of the EOS table, we show the energy and 
compositions as functions of baryon density by choosing T = 10 MeV and Yc = 0.4 
as an example. In the upper-right panel of Fig. [2j we plot energy per baryon (E/B) 
in EOSy together with the results in Shen EOS for comparison. At p B > 0.4 fm~ 3 , 
the energy is lower in EOSV than in Shen EOS due to hyperons (See Fig. [5]). In the 
lower-right panel of Fig. [2} neutron and proton chemical potentials are shown. In EOSY", 
the difference of chemical potentials between neutron and proton is small and neutron 
chemical potential becomes smaller than proton chemical potential at 0.7 fm -3 . This 
is because more SU/(3) symmetric matter is preferred toward high densities, and the 
charge fraction under j3 equilibrium decreases below Yq — 0.4 at high densities. 

We show the particle compositions as functions of baryon density in EOSV and 
EOSY"7r at Yc = 0.4 in the left panels of Fig. [51 In moderately isospin asymmetric 
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Figure 6. Hypcron fraction contours and adiabatic paths in supernova matter at 
Yc = 0.4 from the hyperonic EOS table without pions (EOSF). The contours of the 
fixed number fraction of hyperons (sum of strange baryons) are shown by dashed lines. 
The solid lines denote the contour of fixed entropy per baryon (isentropy) . The dotted 
line shows the trajectory of the dense matter at center during core collapse and bounce. 



matter (Yc = 0.4), the fraction of A particle grows at around p B ~ 0.4 fm -3 , and 
becomes comparable to nucleons at higher densities. The fractions of other strange 
baryons increase slowly and remains small until very high density. In the right panels of 
Fig. [51 we show the particle compositions at Yq = 0.2. With this small charge fraction, 
total hyperon fraction reaches 1 % at 0.25 fm -3 . 

In Fig. [61 we plot the contour map of the fraction of hyperons (sum of strange 
baryons) in the density-temperature plane. In order to have a significant amount of 
hyperons, one needs high density or temperature. In supernova core, the entropy per 
baryon is typically around 1-2 kg, therefore, one needs high densities 0.3-0.4 fm -3 to 
have 1% mixture of hyperons and 0.45 fm -3 for 10%. 

In supernova matter, the effects of hyperons and pions are limited. When isospin 
asymmetry is not high (ex. Yc = Y e = 0.4), total amount of hyperons which appears at 
0.4 fm -3 is around 1% of baryons. The value of Yc remains high due to the neutrino 
trapping during the collapse and bounce [48]. However, after the deleptonization due to 
neutrino emission, Yc becomes smaller and hyperons may appear in the proto-neutron 
star cooling process. Dense matter at higher densities and temperatures may appear 
also in black hole formations, and hyperon effects can be expected in such processes. 
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3.3. Applications to core-collapse supernovae 

As an application of the EOS table with hyperons, we perform the numerical calculations 
of hydrodynamics of core-collapse supernovae. This calculation is aimed to test the data 
of EOS table for numerical simulations and to provide the basic information on the 
properties of EOS in supernovae such as the appearance of hyperons. For this purpose, 
we calculate the adiabatic collapse of iron core of massive stars of 15M [37]. In the 
same way as the hydrodynamical calculations in Sumiyoshi et al. [TO], we calculate 
the general relativistic hydrodynamics under the spherical symmetry without neutrino- 
transfer, which is time-consuming, by assuming that the electron fraction is fixed to 
the initial value in the stellar model. Numerical simulations by neutrino-radiation 
hydrodynamics are in progress. 

We have found that the adiabatic collapse of 15M star with EOSF leads to a 
prompt explosion. This model explosion is caused by the large electron fraction assumed 
and is quite similar to the case obtained with Shen EOS [10]. The explosion energy is 
almost the same as the case with Shen EOS and the difference turns out to be small 
within 0.5%. We plot the trajectory of density and temperature of the central grid in 
the hydrodynamical calculation in Fig. [6j 

We have examined the appearance of hyperons during the evolution of core-collapse 
and bounce. We find that the fraction of hyperons turn out to be very small within 10 -3 . 
This is because the density does not increase drastically even at the core bounce in the 
current model. The peak density is 0.24 fm -3 which is lower than the threshold density 
0.60 fm~ 3 at temperature 21.5 MeV and electron fraction 0.42, where A hyperons appear 
by the same order as nucleons. This small mixture does not affect largely the dynamics 
in the model explosion. 

A large electron fraction leads to a large proton fraction, and, therefore, suppresses 
the appearance of hyperons [48]. We note here that this is the outcome of simple 
adiabatic hydrodynamics without the treatment of neutrinos. When electron captures 
and neutrino trapping are taken into account [IT] , electron fraction might be smaller 
than the current value and may enhance the hyperon appearance. The hyperons will 
definitely appear in the thermal evolution of proto-neutron stars after 20 seconds [19J, 
during which the central density becomes high and the electron fraction gets smaller. In 
recent findings of black hole formation from massive stars of 40M Q (T2J [13], the hyperon 
EOS is necessary since the density becomes extremely high during the collapse toward 
the black hole. It would be interesting to perform the full simulations of core-collapse 
supernovae and related astrophysical phenomena. 

4. Summary and discussion 

In this paper, we have presented several sets of equation of state (EOS) of supernova 
matter (finite temperature nuclear matter with lepton mixture) including hyperons 
(EOSy) using an SU/(3) extended relativistic mean field (RMF) model with a wide 
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coverage of density, temperature, and charge fraction. Supernova matter EOS is one of 
the most essential parts in numerical simulations of core collapse supernovae. At present, 
two sets of supernova matter EOS (Lattimer-Swesty EOS [15] and Shen EOS p3]) are 
widely used. The constituents in these EOSs are nucleons and nuclei, then it is desired 
to include hyperons, which are believed to appear at high densities. Here we have 
extended the relativistic EOS by Shen et al. [16] by introducing hyperons. 

We start from the RMF parameter set TM1 for nucleon sector [17], which well 
describes the bulk properties of nuclei in the wide mass and isospin range. For hyperon- 
meson coupling constants, we adopt the values in Ref. [2] as the starting points. 
Hyperon-vector meson couplings are fixed based on the flavor-spin SU(6) symmetry, 
and hyperon-scalar meson couplings are determined to give the hyperon potentials 
in nucleonic and hyperonic matter. Hyperon potentials in nuclear matter around the 
normal density, U^\po), are accessible in hypernuclear production reactions. Recent 
developments in hypernuclear physics suggest the following potentials for £ [21], [22], [23] 
and S baryons [21 [25] [26] [27] 

U§P(po) ^ +30 MeV , U^ N) (p ) ~ -15 MeV . (11) 

These potentials are consistent with those in the quark-cluster model for YN 
interaction [37] and a chiral model prediction [39]. In this paper, we have modified 
g a s and g a s to explain these potentials, while other coupling constants are unchanged 
from those in Ref. [2]. 

The £ potential in nuclear matter still has ambiguities. Recent theoretical 
analysis [23] has shown that the shape and absolute values of quasi-free £ production 
spectra are well explained in a Woods-Saxon potential with U-e(pq) ~ +15 MeV. On the 
other hand, a few MeV attractive pocket is known to be required to explain the energy 
shift of £~ atom [3^ 1 [35l 136] . then the central repulsion would be stronger to cancel the 
effects of this pocket. Other theoretical analyses [21] [22] suggest that Us(po) — +30 MeV 
would be preferred in order to explain the shape or the absolute yield in £ production 
spectra. In any of these analyses, £ potential should be repulsive or less attractive than 
that for A, then the effects of £ hyperons are much smaller than those in the attractive 
case. It is to be noted that the ambiguities in Us do not affect the supernova EOS very 
much as far as £ hyperon fraction is small. 

Formation of finite nuclei at low densities is another important ingredient in 
supernova simulations. In the present EOS, effects of finite nuclear formation are 
included by using the Shen EOS [16], in which formation of finite nuclei is included 
in the Thomas-Fermi approximation. Effects from finite nuclei are evaluated by the 
difference of free energy and its derivatives in the Shen EOS from the EOS of uniform 
nucleonic matter (TM1) without hyperons at each (T, p B ,Yc). 

We have examined the properties of the EOS with hyperons in neutron star matter 
(T = 0, /^-equilibrium) and supernova matter. Hyperon effects are significant in neutron 
stars as discussed already in the literature [H [2] [3] H] [5] [6] [7J . Hyperons appear at around 
p B ~ 2p in cold matter under /3-equilibrium and soften the EOS. The maximum mass 
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of neutron stars decreases from 2.17M to 1.55M and 1.63M Q when hyperons are 
included with attractive and repulsive hyperon potentials, respectively. In prompt phase 
in supernova explosions, on the other hand, hyperon effects are found to be small in 
a spherical, adiabatic collapse of a 15M star by the hydrodynamics without neutrino 
transfer. In the case with Yq = Y e = 0.4 as a typical example, hyperon fraction becomes 
meaningful (Yy > 1%) at p B > 0.4 fm~ 3 or T > 40 MeV. In the spherical and adiabatic 
core collapse calculation of a massive star with the 15M (47], the maximum density 
and temperature are found to be (p s , T) = (0.24 fm -3 , 22 MeV), which do not reach the 
region of the above hyperon mixture region. It should be noted that this conclusion is 
model dependent. Hyperons may appear more abundantly in more realistic calculations 
with neutrino transfer, which are in progress. 

We have also discussed the roles of pions in neutron stars and supernovae. In this 
work, we have examined the effects of free thermal pions p]. In neutron star matter, the 
absolute value of the charge chemical potential He = H P — Hn is calculated to be larger 
than the pion mass at p B > p , thus charged pions can condensate as far as the pion- 
nucleon interaction is not very repulsive. The EOS softening from pions is moderate and 
limited in the density range p B < 0.88 fm~ 3 without (p-wave) irN attraction, then the 
maximum mass of neutron stars (1.65M ) is almost the same as that without pions. In 
supernova explosions, temperatures are not very high and pion contributions are small. 
At higher temperatures as in the case of black hole formation or high energy heavy-ion 
collisions, the role of pions should be significant. 

There are several points to be improved for deeper understanding of supernova 
matter EOS. First, it is necessary to examine the coupling constants of hyperons with 
hidden strangeness mesons, ( and <fr, which critically decide YY interaction. In this 
paper, we have adopted g^y and g^y in Ref. j2], where the couplings are determined 
based on the SU(6) relation and a conjecture on the hyperon potential depth in hyperon 
matter. An alternative way to determine these couplings would be to invoke various 
hypernuclear and hyperon atom data, such as the double A hypernuclear binding energy 
in aaHc [49] and atomic energy shifts in E~ atom [341 [35J [36] . In Ref. [36], Tsubakihara 
et al. have determined the scalar couplings of g a \, g^\, g a ^ and g^, by using the double 
A hypernuclear bond energy and the atomic energy shift of S~ atom, while vector 
couplings are fixed from the SU/(3) relations. At present, available data are so scarce 
that we cannot fix these couplings based on the data unambiguously, but future coming 
J- PARC and FAIR facilities will provide much more data on YY interaction. Next, it 
is desired to respect chiral symmetry in order to describe very dense matter, in which 
spontaneous broken chiral symmetry will be partially restored. A chiral symmetric 
RMF model [50] is recently developed based on a scalar meson self-energy derived in 
the strong coupling limit of lattice QCD, and it describes binding energies and radii of 
normal nuclei in a comparable precision to TM1. An SU/(3) extended version of this 
chiral RMF is now being developed [36] • Finally, distribution of finite nuclear species 
may be important at low densities [5T] . At finite temperatures, the entropy increase 
by the formation of various fragments will contribute to gain the free energy compared 
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with the single heavy-nuclear configuration assumed in the Thomas- Fermi approach. It 
is not straightforward but challenging to include nuclear statistical equilibrium (NSE) 
distribution in a consistent way in the EOS based on RMF. 

These challenging developments of hadronic and nuclear physics are important to 
understand the extreme conditions in compact objects and to clarify the mechanism of 
explosive phenomena in astrophysics. 
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Appendix A. Note on the EOS table 

Appendix A.l. Locations of data tables 
The data tables are available on 



\protect\vrulewidth0p1http://nucl.sci.hokudai.ac.jp~chikako/EOS/index.html 



or upon request to A. Ohnishi. On the web page, the tables under the name of 
'***.tbl' (EOS table) and '***.rat' (composition table) are available for the sets using 
= —30, 0, +30, +90 [MeV] at normal density. The most recommended potential is 
= +30 MeV. For other hyperons, we adopt the potential depth Uj^ = —30 MeV 
and lli N ^ = — 15 MeV as described in section 12+1 The EOS table with thermal pions are 
available in addition to the standard choice without pions. As described in section 12.21 
the set with thermal free pions are aimed only for the assessment of pion contributions 
in a simple treatment. Further careful treatment of pion interactions is necessary. 



Appendix A. 2. Definition of quantities in the EOS table 

We list the definitions of the physical quantities tabulated in the EOS table, '***.tbl'. 
We note that the order of quantities in the list is partly different from the original Shen 
EOS table [16] since the current table contains lepton and photon contributions. The 
definitions of quantities follows the ones in Shen EOS unless stated specifically below. 

(1) Logarithm of baryon mass density: log 10 (p B ) [g/cm 3 ] 
The baryon mass density p B is defined by 

p B = M u n B (A.l) 
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where M u and Ub are the atomic mass unit and the baryon number density, 
respectively 

(2) Charge ratio: Y C C = ric/n B The charge density ric is defined by 

no = y^qiUi , (A. 2) 

B 

where qi and rij are the charge and the number density of the baryons and the sum 
runs over the baryon octet. 

(3) Entropy per baryon: S/B [ks] 

The entropy per baryon contains the contributions from hadrons, leptons and 
photons. 

(4) Temperature: T [MeV] 

(5) Pressure: P [MeV/fm 3 ] 

The pressure contains the contributions from hadrons, leptons and photons. 

(6) Chemical potential of neutron: fi n [MeV] 

The chemical potential of neutron is measured relative to the nucleon mass 
Myv = 938 MeV. It is connected with the baryon chemical potential as 

Hn = Hb ~ M N (A. 3) 

(7) Chemical potential of proton: fi p [MeV] 

The chemical potential of proton is measured relative to the nucleon mass Mjv- 
The relation to the baryon and charge chemical potentials reads 

Hp = Ub + A*c - M N (A.4) 

(8) Chemical potential of electron: \i e [MeV] 

The chemical potential of electron is determined by the charge neutrality as 

n e = ^qiUi. (A. 5) 

B 

See below for the descriptions on the lepton contributions. 

(9) Free neutron fraction: Y n 

In the uniform matter, the free neutron fraction is simply the ratio, n n /riB- For 
the non-uniform matter at low density, the definition follows the one in Shen EOS. 

(10) Free proton fraction: Y p 

The ratio, n p /riB, as in Y n above. For the fraction of strangeness baryons, see the 
description below on the table of number fractions. 

(11) Mass number of heavy nucleus: A 

The values of from (11) to (14) are taken from Shen EOS table or zero set above 
normal nuclear density. 

(12) Charge number of heavy nucleus: Z 

(13) Heavy nucleus fraction: Xa 
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(14) Alpha-particle fraction: X a 

(15) Energy per baryon: E/B [MeV] 

The energy per baryon is defined with respect to the free nucleon mass and 
contains the contributions from hadrons, leptons and photons. 

(16) Free energy per baryon: E/B [MeV] 

The free energy per baryon is defined with respect to the atomic mass unit M u and 
contains the contributions from hadrons, leptons and photons. 

(17) Effective mass: M* [MeV] 

The effective mass of nucleon is obtained in the RMF theory for uniform matter. 
In non-uniform matter, we replace the effective mass by the free nucleon mass 
M N . 

Appendix A. 3. Data table of composition 

In order to provide the information on the appearance of hyperons other than nucleons, 
we prepare a separate data table (***.rat) for the number fractions. The number 
fraction, Vj = Ui/nB, is given as a function of (p B ,T,Y c ) in the following order. 

(1) Logarithm of baryon mass density: log 10 p s [g/cm 3 ] 

(2) Temperature: T [MeV] 

(3) Charge ratio: Yq 

(4) Neutron ratio (including neutrons in alpha): Y n 

(5) Proton ratio (including protons in alpha): Y p 

(6) A ratio: Y& 

(7) S~ ratio: Y%- 

(8) S° ratio: V E o 

(9) £+ ratio: V E+ 

(10) H~ ratio: Y s - 

(11) E° ratio: Y E o 

Appendix A. 4- Treatment of leptons 

We describe briefly on the contribution of leptons (electrons, muons and neutrinos) in 
the current study. 

We remark that we take into account muons in the case of cold neutron stars. Muons 
appear abundantly at densities higher than the threshold when the electron chemical 
potential exceeds the muon rest mass. The chemical potentials of electrons and muons 
are related with 



(A.6) 
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This is because neutrinos freely escape from the neutron star and are not trapped 
inside. Accordingly, the contributions of electrons and muons are taken into account in 
the discussions of cold neutron stars. 

For the EOS table for core-collapse supernovae, we add the contributions of 
electrons, positrons and photons while muons are not added because of the following 
reason. In the supernova cores, neutrinos are trapped inside the supernova core and 
the Fermi energies of neutrinos becomes high and non-zero. If the chemical equilibrium 
holds, the chemical potentials follows the relations, 

/V - /^ =He~ Vu e , (A.7) 

n a + = . (A.8) 

The latter relation comes from the fact that the net //-type lepton number is zero. 
Because of positive values of //„ e , r.h.s of the chemical equilibrium is reduced. In 
addition, the appearance of muon requires the production of anti-neutrinos of //-type and 
leads to the negative value of [/,„ . Therefore, the appearance of muons are suppressed 
in supernova core. 

Contributions of neutrinos are not added because the treatment depends on the 
density region. In the central part of supernova core, neutrinos are trapped and 
the chemical equilibrium are reached together with neutrinos. Outside the neutrino 
trapping surface, typically at ~ 10 11 g/cm 3 , neutrinos escape freely and neutrinos do 
not contribute. They also depend on the method of neutrino-radiation in numerical 
simulations. 
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